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The total energy distributions (TEDs) of the emission currents in field emission and surface photofield 
emission and the overlayer-induced modifications in the surface electronic structures fi^om the technolo- 
gically important W surfaces with the commensurate W(100)/Na c(2x2), W(110)/Na (2x2) and 
W(l 1 1)/Na (1x1) overlayers are calculated. The TEDs obtained by our recent numerical method that ex- 
tends the full-potential linear augmented plane wave method for the electronic structures to the study of 
field and photofield emission are used to interpret the shifts of the peaks in the experimental TEDs in 
field emission and photofield emission from the W(IOO) and W(l 10) surfaces at sub-monolayer and mo- 
nolayer Na coverage. Hybridization of the 3s Na states with the pairs of (i/-like surface states of the 
strong Swanson hump in clean W(IOO) and surface resonances in clean W(l 1 1) below the Fermi energy 
shifts these W states by about -1.2 eV and -1.0 eV, thus stabilizing these states, to yield new strong 
peaks in the TEDs in field emission and photofield emission ixom W(100)/Na c(2x2) and W(lll)/Na 
(1x1) respectively. The effect of Na intralayer interactions are discussed and are shown to shift the 
strong s- and/'-like peaks in the surface density of states of W(l 10) below and above the Fermi energy 
respectively to lower energy with increased Na coverage, in agreement with experiments. 



I. INTRODUCTION 

Field emission has been extensively used in 
technological applications for imaging as well as to study 
the electronic properties at metal and semiconductor 
surfaces and thin film overlayers''^. Modifications in the 
electron emission characteristics of transition metals due to 
the adsorption of alkali metal overlayers and their oxides is 
of great interest both fundamentally and technologically^. 
The reduction in the work function, induced by the 
overlayer, enhances the quantum mechanical tunneling of 
electrons through the surface making these interfaces 
suitable for thermionic cathode and photoelectron emission 
device applications. To understand and improve the 
performance of these devices, it is important to 
microscopically interpret the physical mechanisms at the 
interfaces that modify the electronic structures and 
emission current with changing coverage. Details of the 
theoretically extracted electronic structure are expected to 
shed light on the role played by surface states and surface 
resonances in the substrate - overlayer bonding, the nature 
of bonds formed, the charge shifts and modification of the 
potential induced at the interface by the overlayer. While 
several such overlayer systems have been extensively 
studied theoretically"''^'^, the theoretical work on the 
tungsten-sodium (W-Na) interface has been limited to only 
the (100) plane^, or was using simplified formalisms that 
neglect the dependence of the substrate electronic structure 
on surface orientation^'**. 
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Most experimental studies of the W-Na interfaces 
were focused on work function measurements and under- 
standing the Na growth mechanism^' ^"'^. Photoemission 
measurements'^ from W(llO) covered by Na, K and Cs 
from low up to complete one monolayer coverage were 
used to analyze the binding energies of the alkali s, and 
tungsten 3p and 4f core electrons. Metastable Impact Elec- 
tron Spectroscopy and Ultraviolet Photoemission Spectros- 
copy were used to study the electron spectra from Li, Na, K 
and Cs overlayers on W(llO)'''. An experimental investiga- 
tion of the valence electronic structure of the W-Na- 
vacuum interfaces for W(IOO) and W(llO) was carried out 
by Derraa and Lee'^ by the surface sensitive field emission 
(FE) and photofield emission (PFE) techniques. Na cover- 
age was in the range from to 9= 1/2 and to 9= 1/4 re- 
spectively, where 9 is the ratio between the number of 
atoms in the overlayer and that in the outermost substrate 
layer. 

Although many experimental results on W-Na have 
been available for more than a decade and in spite of the 
technological importance of these interfaces, only a few 
theoretical investigations - including only one recent^ pa- 
rameter-free, ab-initio study- of the electronic structure of 
these interfaces are reported. The recent work^ that investi- 
gated Na overlayer growth at W electrodes at different Na 
pressure in high pressure Na discharge lamps discusses the 
pressure dependent work function at the W(IOO) surface at 
different coverages. The measured work function was com- 
pared with results extracted from electronic structure calcu- 



lations of W(IOO) with commensurate Na overlayers at dif- 
ferent coverage using the pseudopotential plane wave me- 
thod based on density functional theory (DFT). Lang^ car- 
ried out surface electronic structure calculations of a metal 
surface with a single Na atom, and Ishida^ calculated the 
electronic structure of Na interfaces at different Na cover- 
age, both applying the jellium model to the substrate. These 
results^' ** were used in Ref. [15] to interpret the experimen- 
tal FE and PFE data. The low-index surfaces of W, howev- 
er, contain several d-like surface states and surface reson- 
ances'^''^ that are crucial for the understanding of the sub- 
strate -overlayer bonding, and which the jellium model was 
not able to accurately describe. 

In the present paper, we report the modifications in the 
surface electronic structures of W(IOO), W(llO) and 
W(lll) substrates due to commensurate Na overlayers 
using the ab-initio full-potential linear augmented plane 
wave (FP-LAPW) method in order to interpret the existing 
experimental FE and PFE data'^. FE and PFE are very 
sensitive to the electronic states close to the center of the 
surface Brillouin zone (SBZ), because the surface potential 
barrier attenuates exponentially the electronic states of 
increasing wavevector. To accurately compare the 
theoretical results with the experiments we extend the 
electronic structure calculations to yield the total energy 
distributions of the emission currents (TEDs) in FE and 
surface PFE and the A:-resolved layer densities of states 
(K-LDOS). 

Low-energy electron diffraction (LEED) 

measurements" demonstrate that Na from an atomic beam 
grows pseudomorphically layer by layer on W(100)/Na up 
to 80 layers in a c(2x2) structure corresponding to ^= 1/2. 
Our W(IOO) calculations were carried out with a Na c(2x2) 
overlayer [denoted W(100)/Na c(2x2)], which corresponds 
to one monolayer (ML) coverage. On W(l 10), LEED data'" 
show that the Na overlayer undergoes a number of 
successive commensurate structures before completing the 
first atomic ML, which is an incommensurate, hexagonal 
layer at a coverage corresponding to ^ = 3/5. Since our 
numerical method, based on the repeated supercell 
geometry, involves the assumption of translational 
symmetry parallel to the surface, only commensurate 
overlayers were considered. We carried out emission 
current calculations at an commensurate overlayer 
corresponding to 6* = 1/4 (denoted W(110)/Na (81,4) , S 
stands for structure), which represents a (2x2) overlayer 
observed in the LEED data'", in order to compare them 
with available experimental results'^. To study the 
modification of the electronic structure with increased Na 
concentration we also report results of W(llO) substrates 
with overlayers at 6*= 2/5 [W(110)/Na (S2/5)], which is the 
commensurate overlayer of highest density observed by the 
LEED measurements, and at 9 =1/6 [W(110)/Na (Sye)], 
which is the lowest-coverage commensurate overlayer 
observed. Regarding the Na overlayer on W(lll), despite 



that W has a smaller atomic radius than Na, the low 
atomic density of the W(lll) substrate allows it, in 
principle, to accommodate a Na (1x1) overlayer of equal 
atomic density. Our W(lll) calculations were carried out 
with a hexagonal Na (1x1) overlayer [denoted W(lll)/Na 
(1x1)] To the best of our knowledge, no experimental data 
on the structure of Na overlayers on W(l 11) exists. 

The electronic structure calculations are outlined in 
Section 11. The method of calculation has been described in 
more detail elsewhere'^'. In Section III the calculations of 
W(100)/Na c(2x2) are interpreted and compared to the 
experimentally measured TEDs in FE and PFE. Similar 
discussions of W(110)/Na (52/5), W(110)/Na (^1/4), 
W(110)/Na (5i/6) and W(lll)/Na (1x1) follow in Sections 
IV and V respectively. The results and conclusions of this 
work are summarized in Section VI. 

II, COMPUTATIONAL METHODS 

The metal-adsorbate-vacuum interface is represented 
by a supercell, whose electronic structure is calculated self- 
consistently on the basis of density fiinctional theory by the 
FP-LAPW method. Exchange and correlation are treated in 
the generalized-gradient approximation'**, and relativistic 
corrections including the spin-orbit interaction are taken 
into account unless otherwise specified. We used the 
electronic structure software package WIEN2K'^, which 
was modified as described in Refs. 16 and 17 to calculate 
the TEDs in FE and surface PFE, and the ratio of the TED 
of the metal to that of the free-electron metal (the 
enhancement factor of the emission current). We also 
report the A-resolved surface density of states (K-SDOS), 
which describes the surface density of states (SDOS) 
weighted by the normalized tunneling factor, hence 
emphasizing the features of the electronic structure in the 
vicinity of the center of the SBZ (F). 

We found that for satisfactory convergence of our sur- 
face electronic structure calculations of clean W, the super- 
cell must have at least two intermediate layers in register 
with the surface and the bulk. The geometry of the 
W(100)/Na c(2x2) interface is therefore described by a su- 
percell that consists of 13 W layers stacked parallel to the 
( 1 00) plane, surrounded on each side by a region of half 
that volume containing one Na atom. The W(l 10)/Na (S1/4) 
supercell, corresponding to a (2x2) overlayer, consists of 
13 W layers stacked parallel to the (110) plane, and the 
W(l 1 1)/Na (1x1) supercell consists of 19 W layers stacked 
parallel to the W(l 1 1) plane surrounded on each side by a 
region of half that volume containing one Na atom. 

In Fig. la, 2a and 3a the solid square and rhombuses 
represent primitive unit cells of the Na c(2x2), (2x2) and 
(1x1) overlayers, and the dashed square and rhombuses 
represent primitive unit cells of the clean W(IOO), W(llO) 
and W(lll) substrates respectively. The corresponding 
SBZs are shown in Fig. lb, 2b and 3b respectively. The 



symmetry points of the SBZ of the overlayers are denoted 
by a prime. The effect of the c(2x2) overlayer is to fold 
back the symmetry point M of the (1x1) substrate to r' in 
the SBZ of the overlayer. The effect of the (2x2) overlayer 
is to fold back the symmetry points N and S of the sub- 
strate to r' in the SBZ of the overlayer. 



underestimate the strengths of some peaks in the TEDs in 
PFE that are predominantly due to bulk photoexcitation. 
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Fig. 1. (Color online) Surface unit cells and surface Brillouin 
zones of clean W(IOO) and of a Na c(2x2) overlayer on W(IOO). 

(a) Atomic positions in the top two layers of W(IOO) and in the 
overlayer. Also shown are the primitive unit cells of the clean 
substrate (blue dashed lines) and of the overlayer (red solid lines). 

(b) The surface Brillouin zones of the substrate and the overlayer 
plotted in the correct orientation relative to the unit cells in (a). 

In interpreting the PFE results, it is important to 
distinguish between electron transitions resulting from p 
and s polarization of the incident light. In p polarization, 
electrons are excited by surface photoexcitation, so the 
final states are a continuum of ixee-electron-like states just 
outside the surface. Hence any features in the calculated 
TEDs in PFE in surface photoexcitation correspond to the 
initial states of the optical transition. If the polarized beam 
is also partially s polarized, electrons are also excited by 
bulk photoexcitation, hence direct transitions are induced 
between the initial and final states inside the metal. Hence 
any features observed in bulk photoexcitation correspond to 
the initial and final states of the optical transition. Because 
the available experimental PFE data are carried out 
predominantly in p polarization and because the transition 
rate in surface photoexcitation is typically much stronger 
than that in bulk photoexcitation we limited our PFE 
calculations to surface photoexcitation. This might 
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Fig. 2. (Color online) Surface unit cells and surface Brillouin 
zones of clean W(llO) and of a Na (5i/4) overlayer on W(llO). 
Color and line style descriptions are the same as in Fig. 1 . 

In this paper, the experimentally-observed emission 
peaks are labeled alphabetically and the emitting facet is 
denoted by a numerical subscript (1 denotes 100, 2 denotes 
110 and 3 denotes 111). The calculated peaks are denoted 
by a prime. All energies are expressed relative to the Fermi 
level Ep. 

III. RESULTS AND DISCUSSION FOR W(100)/Na 

III.l. Field and photofield emission currents from 
W(100)/Na c(2x2) 

Experimental TEDs in FE and PFE from clean W(IOO) 
at 70 K^" and at room temperature'^'^'^'^' show a strong 
peak Bi, known as the Swanson hump^^. The calculated 
dispersion plot of clean W(IOO) plotted in the SBZ of a 
c(2x2) structure is shown in Fig. 4a. The calculated K- 
SDOS of clean W(IOO), plotted in Fig. 5a, shows a peak B/ 
that is due to a pair of surface states B/ of (//-like symme- 
try close to F' (Fig. 4a) that are responsible for the ob- 
served strong peak Bi, where z denotes the direction along 
the normal to the surface. The energies of the calculated 
and experimental peaks are compared in Table I, and the 
symmetries of the calculated electron states that dominate 
the emission are reported. A detailed analysis of FE and 
PFE Irom clean W(IOO) can be found in Ref 16. 
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Fig. 3 . (Color online) Surface unit cell and surface Brillouin zone 
of clean W(lll) and of a Na (1x1) overlayer on W(lll). (a) 
Atomic positions in the top two hexagonal layers of W(lll) and 
in the overlayer. Also shown is the primitive surface unit cell of 
the clean substrate and of the overlayer. (b) The surface Brillouin 
zone plotted in correct orientation with respect to the correspond- 
ing unit cell in (a). 

The logarithms of the calculated and the experimen- 
tal'^ TEDs of the PFE current from W(100)/Na c(2x2) at 
room temperature as a function of the initial state energy 
are shown in Figs. 5c and 5d (upper curve at 1 ML cover- 
age) respectively. In our calculations the energy of the pho- 
tons in a /^-polarized beam was 3.05 eV and the electric 
field strength was 0.16 V.A"' as determined from the expe- 
rimental data. The energies and symmetries of the peaks in 
FE and PFE are reported in Table 1 . The band structure of 
W(100)/Na c(2x2) along F'X'are shown in Fig. 4c and the 
band structure of an isolated Na layer (Fig. 4b) calculated 
using a supercell in which Na atoms occupy the same sites 
as in the W(100)/Na c(2x2) supercell and the W sites are 
empty is shown in Fig. 4b. When a Na c(2x2) overlayer is 
adsorbed on the clean W(IOO) substrate the s-like states of 
the overlayer near F' (Fig. 4b) hybridize with the surface 
states Bi' of the substrate at -0.27 eV and -0.30 eV (Fig. 
4a), shifting them to -1.49 eV (Fig. 4c) and yielding peaks 
Fi' in the K-SDOS (Fig. 5b) and in the TED in PFE (Fig. 
5c). This is consistent with the observed suppression of the 
experimental peak Bj in the TED in FE from clean W(IOO) 
above a Na coverage of 0.5 ML, and the strong enhanced 
emission in the TED in FE at 1 ML at -1.5 eV, the lowest 
sampled energy, confirming a strong peak Fi at or below 
-1.5 eV'^. In Fig. 5d, the experimental TED in PFE at 



1 ML shows a strong peak Fi that is partially attributed to 
peak Fi' calculated in surface photoexcitation. 

Clean W(l«0)c(2x2) Isolated Na c(2x2) VV (100) /INa 0(2x2) 








I A K H,. 



/ ^v;-(&J jpj|il^:=";|:;='!i:.j[ci 



..'\ 



•:■■ .■!ii!''.;i^^-':iK.e*' 









\.,.^-/f^i^'^ix^. 




.^''^U' Ji- ';■- — ■■- 



'i^i^:M^^' 



'- K 



K VECTOR K VECTOR K VECTOR 

Fig. 4. (Color online) Dispersion along X F of (a) clean 
W(IOO), (b) an isolated Na c(2x2) overlayer, and (c) W(100)/Na 
c(2x2), all plotted in the SBZ of the overlayer. Surface states and 
surface resonances are shown by thick curves. The predominant 
symmetry in the surface layer is shown by the line style {s, red 
(dashed);/), green (dotted); and d, blue or grey (sohd)]. Bulk and 
intermediate states are shown by thin grey (dotted) curves. 

In Fig. 4a the bands of clean W(IOO) are plotted in the 
SBZ of the c(2x2) overlayer in order to visualize the energy 
shifts when comparing them with the bands of W( 1 00)/Na 
c(2x2). As a result the symmetry point M in the SBZ of 
clean W(IOO) is folded back to F' in the SBZ of the over- 
layer. This effect is responsible for the d^Aiks surface 
states Ci' with energies +1.7 eV and +2.18 eV at F' in Fig. 
4a. Because FE is dominated by emission from electron 
states close to F', these states at M of the (1x1) W(IOO) 
substrate are expected to make a negligible contribution to 
emission from clean W(IOO) but can make a significant 
contribution to emission from W(100)/Na c(2x2) in bulk 
photoexcitation. In W(100)/Na c(2x2), the/'z-like states of 
the overlayer hybridize with the surface states C/ shifting 
them to lower energy to yield the strong peaks J/ in the K- 
SDOS. The experimental enhancement factor in FE at 
1 ML coverage shows a strong peak J i at or above +1.0 
eV'^ that is attributed to the surface states Ji'. The exact 
energy of peak J i cannot be determined because the expo- 
nential cutoff in the TED in FE that is due to the reduced 
electron occupation above Ep prevents analyzing the data 
above +1.0 eV. Because the energies of states F/ and J/ 
differ by roughly the photon energy used, it is expected that 
bulk photoexcitation between these states will contribute to 
the strength of the experimental peak Fi. The strong peak 
Fi observed in PFE is therefore attributed to surface pho- 
toexcitation from the initial states F/, as well as bulk pho- 



toexcitation between the initial states F/ and final states J/. 
Peak Ai' in the calculated K-SDOS of clean W(IOO), due to 
a region of dxz+dyz-like surface resonances Ai' at about 
-0.7 eV close to F', is responsible for the weak peak Ai 
observed in the experimental TED in FE^^. While peak Ai' 
and Bi' are comparable in the K-SDOS, peak Ai' is much 
weaker in the TED in FE than peak B/, consistent with 
experiments, due to the different symmetry of the underly- 
ing electron states. In FE only the term of smallest parallel 
wavevector (ky+Gn) in the plane wave expansion of the 
wave function contributes effectively to the emission cur- 
rent; k[| denotes a wavevector in the first SBZ and Gy de- 
notes a reciprocal lattice vector'^'^. While this term is sig- 
nificant in S-, p^-, and d^Aike states at r , because the wave 
functions of these states are symmetric in a plane parallel to 
the surface, it is small or zero in Px+Pv-, d^z+dy^-, d^. ^- and 
c4j,.-like states, because the wave functions are antisymme- 
tric in a plane parallel to the surface. The 5-like states of the 
overlayer hybridize with the dxz+dyz-like surface resonances 
Ai' in an extended region along F'X' to yield peak Gi' in 
the calculated TED in PFE. The experimental TED in PFE 
at a Na coverage from 0.4 to 0.8 ML'^ show a strong peak 
Gi (see lower curve in Fig. 5d at 0.6 ML coverage) that is 
attributed to the electron states G/. 

Ref. 15 also suggests a peak labeled D at a final state 
energy at +3.2 eV in the TED in PFE that was attributed to 
a high DOS at the final state energy. In Fig. 5 of Ref 15, 
however, the TED in PFE with 2.61 eV photons at 1 ML 
was shifted erroneously in energy by +0.5 eV resulting in 
the suggestion that peak D was due to a high DOS at the 
final state energy, i.e. due to bulk photoexcitation. Our cal- 
culations show that at x' , the px+Pv -like states of the over- 
layer hybridize with the d.^ -like surface states Di' of clean 
W(IOO) above Ep shifting them to -0.22 eV. The exten- 
sions of these bands of surface states and surface reson- 
ances that disperse along X' F' to higher energy crossing Ef 
yield peak H/. The weak peak Hi (labeled D in Ref. 15) in 
the experimental TED in PFE at 1 ML Na coverage is attri- 
buted to peak Hi' calculated in surface photoexcitation. 

III.2. Layer densities of states of W(100)/Na c(2x2) 

The layer densities of states (LDOSs) in the central (bulk) 
layer and in the surface layer (SDOS) of the supercell of 
clean W(IOO), shown in Figs. 6a and 6b respectively, show 
large differences due to the breaking of translational sym- 
metry and the reduced atomic coordination at the surface. 
Strong peaks in the bulk layer are suppressed at the surface, 
and strong surface resonance peaks at about -4.1 eV, 
-0.3 eV and +1.6 eV are suppressed in the bulk. 

The LDOS of an isolated Na layer, calculated using a 
supercell in which Na atoms occupy the same sites as in the 
W(100)/Na c(2x2) supercell and the W sites are empty, is 
shown in Fig. 6c. The LDOS is nearly ixee-electron-like 
below Ef. From -1.2 eV to +1.0 eV it is predominantly 5- 



like with minor admixture of /'.v+fv-like states and above 
+1.0 eV it is predominantly /»-like. The LDOSs in the sub- 
strate and in the overlayer of the supercell of W(100)/Na 
c(2x2) are shown in Figs. 6d and 6e respectively. Due to 
hybridization of the s- and p-like states of the overlayer 
with the states of the substrate the LDOS in the overlayer 
deviates strongly from its nearly free-electron-like nature in 
the isolated state. It spreads out to -4.5 eV and shows 
strong 5-like peaks at -4.1eV and -1.5 eV, and /'v+Pv-like 
peaks at -0.3 eV and +1.1 eV due to the interaction of the 
Na states with the W surface resonances. The strong sur- 
face resonance peaks of W(IOO) below Ep are only slightly 
modified thereby. 
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Fig. 5. (Color online) (a) K-SDOS of clean W(IOO) and (b) K- 
SDOS of a Na c(2x2) overlayer on W(IOO). The successive 
curves in the cumulative plots (a) and (b) show the contributions 
of s-like (red), p-like (green), t/^-j +t//./-like (grey), dxz+dyz-like 
(light blue) and dz -like (dark blue) states, respectively. The shad- 
ing denotes the dz -like contributions, [(c) and (d)] TEDs in PFE 
for W(100)/Na c(2x2) with 3.05 eV photons, plotted as a function 
of the initial state energy. The calculated plot (c) is based on sur- 
face photoexcitation. The experimental plot (d) shows a strong 
peak Fi that is due to surface as well as bulk photoexcitation. 

In summary, our calculated results are consistent with 
the peaks observed in FE and PFE from clean W( 1 00) and 
W(100)/Na c(2x2) (Table I), indicating that the present cal- 
culations give a realistic picture of the changes in the elec- 
tronic structure of W( 1 00) due to adsorption of a Na c(2x2) 
overlayer. 

For the muffin-tin radius of 1.7 A used in the present 
calculation the total valence charge of 1.00 electrons within 
the muffin-tin sphere of Na does not change significantly 



Table I. Comparison of the peaks observed in the experimental TEDs in FE and PFE from clean 
W(IOO) (Ai and BO'"'' and W(100)/Na c(2x2) (Fi to JO'', with peaks in the calculated TEDs. s.s. 
stands for surface states and s.r. for surface resonances. 



Experimental 


Calculated 


Peak Label E-Ep (eV) 


Peak Label E-Ep (eV) Character Symmetry in 

Overlayer [Substrate] 


Ai 
Bi 


-0.73 (5) 
-0.32 (3) 


Ai' 
Bi' 


-0.69 (2) s.r. 
-0.32 (2) 5.5. 


[dh 


Fi 


<-1.5 


Fi' 


-1.50(2) 


s.s. 


5 [^/] 


Gi 


-0.8 (1) 


Gi' 


-0.63, -0.83 (2) 


s.r. 


5 [d^z+dy,] 


Hi 


Close to +0.0 


Hi' 


-0.05+0.04 (2) 


s.r. 


Px+Py [dxz+dy^] 


Ji 


<+1.55 


Ji' 


+1.2 to +1.6 


s.s. -s.r. 


p. [dz'] 
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Fig. 6. (Color online) LDOS of clean W(IOO) in the (a) central 
(bulk) layer and (b) surface layer, (c) LDOS in an isolated Na 
c(2x2) layer. LDOS of W(100)/Na c(2x2) in the (d) W sub- 
strate and (e) Na overlayer. In these cumulative plots, the 
areas between successive curves show the contributions of s- 
like (red), p-\ike (green), t/j^.+cZ/.v^-like (grey), d^^+dyz-Mke 
(light blue) and d^ -like (dark blue) states, respectively. 

upon adsorption. Instead, the total charge is redistributed 
among the various symmetry components. The 5-like 
charge (integrated LDOS) decreases by about 0.43 elec- 
trons per Na atom and the p-and d-like charge increase 
by about 0.3 1 and 0. 12 electrons per Na atom respective- 
ly. Hence, nearly no charge transfer takes place from the 
overlayer to the substrate. The same conclusion was re- 
ported'^ based on experimental photoemission measure- 
ments from W( 110) with Na, K and Cs overlayers at co- 
verages up to one atomic layer that show little if any 



shifts in binding energies (energy with respect to the 
Fermi energy) of the core 3p and 4f W-electrons. 



III.3. Work function of W(100)/Na c(2x2) 

Typically, valence electrons spill out from the sur- 
face of a clean transition metal, producing an inwardly- 
directed surface electric dipole layer that increases the 
work function. Alkali and alkaline-earth overlayers on 
transition metal surfaces produce additional surface elec- 
tric dipole layers whose net effect is to decrease the 
work function. 

The spatial redistribution of the valence electron 
density that occurs when a Na c(2x2) overlayer is ad- 
sorbed on W( 1 00) was calculated by subtracting the sum 
of the valence electron density distribution of the clean 
W substrate and of the isolated Na c(2x2) overlayer from 
the valence electron density distribution of W(100)/Na 
c(2x2). The redistribution of the electron density in a 
{110} plane that passes through the metal surface is 
shown in Fig. 7. Electrons move out from the regions 
surrounding the W and Na atoms and accumulate in a 
layer between the overlayer and the substrate. The net 
effect is an outwardly-directed dipole layer resulting in a 
Na-induced lowering in the potential energy of an elec- 
tron outside the surface with respect to the bulk and thus 
lowering the work function. In the present work, the 
work function of W(100)/Na c(2x2) is 2.1 eV as esti- 
mated by calculating the difference between the Cou- 
lomb potential energy far into the vacuum region of the 
supercell and the Fermi energy. The calculated work 
function agrees well with the value of 2.2 eV observed 
experimentally at 77 K using a field-emission micro- 
scope', and with the value of 2.3 eV calculated by means 
of the pseudopotential plane wave method based on den- 
sity functional theory ^ 
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Fig. 7. (Color online) Calculated redistribution of the valence 
electron density that occurs when a Na c(2x2) overlayer is 
adsorbed on a W(IOO) surface, plotted in a {110} plane that 
intersects the surface at right angles. The dashed rectangles 
represent the planes of the overlayer (labeled Na) and of the 
substrate (labeled W). A positive shift corresponds to an in- 
creased electron density. 

IV. RESULTS AND DISCUSSION FOR W(110)/Na 

IV.l. Field and photofield emission currents from 
W(110)/Na(5i/4) 

The calculated TEDs in FE and PFE from clean 
W(llO)'^ are consistent with the experimental TEDs at 
78K^° and at room temperature'^''^ in that they show lit- 
tle initial state structure. The observed TEDs in FE in the 
range of Na coverage from ^ = to 1/4 show no struc- 
ture, but above that coverage they show several peaks'^ 
In Table II, the energies and symmetries of the calcu- 
lated peaks in the TEDs in PFE at 3.05 eV photons (Fig. 
8c) and the K-SDOS of W(110)/Na {S^^) (Fig. 8b) are 
listed, and the energies of the peaks are compared with 
those observed in the PFE experiments (Fig. 8d). 

To facilitate the interpretation of the charge shifts in 
W(llO) due to the Na {Syi) overlayer the energy bands 
and LDOS (Fig. 9c) of an isolated Na layer calculated 
using a supercell in which Na atoms occupy the same 
sites as in the W(l 10)/Na (5'i/4) supercell and the W sites 
are empty have been calculated. The 5-like valence states 
form a nearly-free-electron band along F' S' . The bot- 
tom of the band is at -0.79 eV. The p^-Viks states of the 
isolated Na (5'i/4) layer form another band at +1.39 eV 
that disperses to higher energy along F' S' while the 
Pv+Pv-like states form a band at +1.06 that disperses to 
higher energy along S' F' . 



Fig. 8. (Color online) (a) Dispersion plots along S' F' in the 
surface Brillouin zone of W(110)/Na (5i/4). Surface states and 
surface resonances are shown by thick curves. The predomi- 
nant symmetry in the surface layer is shown by the line style 
[s, red (dashed); p, green (dotted); and d, blue (sohd)]. Bulk 
and intermediate states are shown by thin grey dotted curves, 
(b) K-SDOS of W(110)/Na {Sw). The successive curves in the 
cumulative plots show the contributions of s-\\ke (red), Px+Pv- 
hke (dark green), p^-\\ke (hght green) and d-\\ke, (blue) states, 
respectively, [(c) and (d)] TEDs in PFE for W(110)/Na (^1/4) 
with 3.05 eV photons, plotted as a function of the initial state 
energy. The calculated plot (c) is based on surface photoexci- 
tation. The experimental plot (d) shows two additional final 
state peaks K2 and J2 that are attributed to bulk photoexcita- 
tion. 

The TEDs in PFE with 3.05 eV photons in the range of 
coverage irom 6= 1/7 to 6= 1/4 (Fig. 8d) show a strong 
peak K2 that shifts slightly to lower energy with in- 
creased coverage, and that has been attributed to bulk 
photoexcitation'^. At the corresponding final state ener- 
gy of +2.45 eV the K-SDOS of W(110)/Na {S^) shown 
in Fig. 8b shows a strong peak K2'. In the presence of a 
(2x2) overlayer, the symmetry point s in the SBZ of 
clean W(llO) is folded back to F' in the SBZ of the 
overlayer. Bands of px-^py-V\ks, states of the isolated Na 
layer close to F' hybridize with J^-like W states at 
+2.74 eV in clean W(llO) that have been folded back 
irom S to yield the bands of surface resonance states K2' 
(Fig. 8a). The calculated TED shows a weak peak G2' at 
the corresponding initial state energy that is due to a 
band of 5-like surface resonances. The observed peak K2 
is attributed to bulk photoexcitation between these states. 
The observed TED in PFE shows another peak J2 at 
1 ML. Peak J2' in the K-SDOS is due to s and p,-\\ks 
states of the overlayer close to F' that hybridize with d^- 
like states of the substrate at +1.89 eV. The observed 
peak J2 is attributed to bulk photoexcitation to these 
states. 



Table II. Comparison between the energies of the peaks observed in the TEDs in PFE from 
W(l 10)/Na (^im)'^ and the energies of the calculated peaks, s.r. stands for surface resonances. 



Experimental 


Calculated 


Peak Label E-Ep (eV) 


Peak Label 


E-Ef{qW) 


Character 


Symmetry in 
Overlayer [Substrate] 






F2' 


-2.08 (2) 


s.r. 


s [d^,+dy,] 


G2 


-0.60 (3) 


G2' 


-0.60 (2) 


s.r. 


s [d^,+dy,] 






H2' 


+0.80 (2) 


s.r. 


PzIPz] 


h 


+1.65 (3) 


J2' 


+L80(5) 


s.r. 


s.pdd.'] 


K2 +2.45 (3) 


K2' 


+2.5(1) 


s.r. 


P.+Py \d^\ 



The 5-like states of the Na (5i/4) overlayer close to 
F' hybridize with the W states to form bands of 5-like 
surface resonances F2' that yield a strong peak F2' in the 
calculated TED in PFE at -2.08 eV. Because the energy 
of the peak of the surface potential barrier is about 
+1.55 eV in the PFE experiments of Fig. 8d, only elec- 
trons having normal energy greater than -1.50 eV can 
pass above the peak of the barrier which accounts for the 
cutoff below about -1.5 eV. Peak F2' could not be veri- 
fied by the experimental results because its energy range 
is below the cutoff energy in the experimental TED. 

IV.2. Layer densities of states of W(110)/Na 

Clean W(llO) shows prominent peaks in the bulk 
LDOS (Fig. 9a) that are greatly attenuated in the surface 
layer (Fig. 9b) and shows prominent surface peaks at 
-3.0 eV, -1.5 eV and +2.3 eV, that grow progressively 
weaker in the subsurface layers. The LDOS of an iso- 
lated Na layer, calculated using a supercell in which Na 
atoms occupy the same sites as in the W(110)/Na (51/4) 
supercell and the W sites are empty, is shown in Fig. 9c. 
The LDOS is 5-like below +0.8 eV and is well separated 
from the/'-like LDOS that predominates above +1.0 eV. 
In W(l 10)/Na (Sys) overlap between the s and/»-like Na 
states and the W states causes the energy gap in the over- 
layer above Ep (Fig. 9e) to disappear. The LDOS in the 
overlayer spreads out in energy to about -3.5 eV and 
shows a strong 5-like peak at -2.1 eV and several pre- 
dominantly />-like peaks above Ef. The strong surface 
resonance peaks of clean W( 110) shift slightly in energy 
(Fig. 9d) and extend into the overlayer. 

Peak K2 in the experimental TEDs in PFE from 
W(110)/Na shifts to lower energy with increased cover- 
age'^ Electronic structure calculations of a metal sub- 
strate represented by a jellium model with Na overlayers 
arranged in a square lattice at different coverages by 
Ishida**, on the other hand, show no shift in the SDOS to 
lower energy of the valence states with increased Na 
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coverage. To study the effect of Na coverage on the 
SDOS, we report the results of two more W(l 10) calcu- 
lations at different coverages; W(110)/Na (5i/6) at ^ = 
1/6 has the lowest-coverage commensurate overlayer 
observed by LEED data, and W(l 10)/Na (^2/5) at 6*= 2/5 
has the highest-coverage commensurate overlayer ob- 
served'". The structures of the overlayers correspond to 
structures b and fin Fig.4 of Ref. 10 respectively. The 
same normal distances between the Na overlayers and 
the W substrate in the three calculations was used. Due 
to the computational challenges introduced by the large 
number of 8 1 inequivalent atoms in the supercell, spin- 
orbit interaction has not been included in the calculations 
at 9= 1/6. Spin-orbit interaction in clean W(l 10) is 



Clean W(110) 

Center Loyer Surface 




Isolated WC110)/Na (S,y^) 

Na CSi/J W Substrate Ma Overlayer 




\\ 
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Fig. 9. (Color online) LDOS of clean W(l 10) in the (a) central 
(bulk) layer and (b) surface layer, (c) LDOS in an isolated Na 
(5i/4) layer. LDOS of W(110)/Na (^1/4) in the (d) W substrate 
and (e) Na overlayer. In these cumulative plots, the areas be- 
tween successive curves show the contributions of 5-like (red), 
p-\\ka (green), d^y+d/./-like (grey), d^^+dyA'Afs (light blue) 
and d^ -like (dark blue) states, respectively. 



responsible for minor shifts in energy of some energy 
states (a maximum of 0.2 eV for states below Ef), and in 
the much lighter Na they show negligible effect. While 
inclusion of spin-orbit interaction in W(l 10)/Na (5'i/6) is 
expected to cause splitting of some states, the resulting 
energy shifts, as judged by the clean W(llO) case, are 
expected to be smaller than the energy shifts caused by 
the coverage variations discussed in the present work, 
and hence they do not invalidate the discussion and con- 
clusions in the coming paragraphs. 

The SDOSs in the three Na overlayers on W(llO) 
of coverages 0= 1/6, 1/4, 2/5, are shown in Figs. 10a, b 
and c. Overlap between the Na states and the W states 
results in a noticeable p^, py and p^-like component in the 
SDOSs below Ep. Increasing the Na coverage broadens 
the peaks in the SDOS, especially above Ef. It also shifts 
the peaks to lower energy, for example, the s-like peak 
below Ef at -2.0 eV shifts by about -0.4 eV when the 
coverage is increased from 9= 1/6 to 6=115, and the 
strong s and/'z-like peak at +1.7 shifts by about -1.2 eV. 
Also the /'v+Py-like peak at +2.5 eV shift by about 
-0.3 eV when the coverage is increased from 0= 1/4 to 
9= 215, which is consistent with the experimental TEDs 
in PFE from W(110)/Na'^ that show a shift to lower 
energy of peak K2 with increased coverage. 

To separate the effect of intralayer Na-Na interac- 
tions from the interlayer W-Na interactions the LDOSs 
is reported for isolated Na layers using supercells in 
which Na atoms occupy the same sites as in the corres- 
ponding W(110)/Na supercells at coverages 6= 1/6, 1/4, 
2/5, and the W sites are empty (Figs. Ha, b and c re- 
spectively). At low Na concentration intralayer interac- 
tions between the 2>s and 2>p atomic electron states of Na 
on adjacent atoms broaden them into bands. This results 
in a 5-like LDOS extending over a broad energy range 
through Ep and a mainly j3-like LDOS at higher energy. 
While in the isolated Na {Sye) and Na (^i^) layers the s- 
like states are well separated from the /'-like states, in the 
isolated Na (52/5) layer overlap between the atomic 5- 
and /'-like states results in a continuous LDOS, and a 
weak px- and py-Vike component in the LDOS below Ef. 
It is seen that intralayer interactions shift, for example, 
the lowest energy 5-like valence states in the LDOS of 
the isolated Na layers, which are states at F', by -1.0 eV 
when the Na concentration is increased from 9= 1/6 to 
9= 2/5, and explains the decrease in energy of the strong 
5-like peak in the SDOS below Ep with increasing 
coverage. 



W(110)/Na(S,/(,) W(110)/No(S,/^ W(1 1 0)/Na (Sj/5) 




0,5 10 0.5 10 0.5 

LDOS (Electrons Atom"^ eV"') 

Fig. 10. (Color online) LDOS of W(110)/Na in the Na over- 
layer at increasing Na coverage for (a) W(110)/Na (Si/e), 
(b) W(110)/Na (Sw) and (c) W(110)/Na {82/5). In these cumu- 
lative plots, the areas between successive curves show the 
contributions of j-Iike (red), p^+py-\\ke (dark green), /»^-Iike 
(light green) and cZ-like (blue) states, respectively. 



No (S^b) 



Isolated Na Layer 
No (S,/^) 




024602460246 
LDOS (Electrons Atom"' eV"^) 

Fig. IL (Color online) LDOS in isolated Na layers at increas- 
ing concentration for (a) Na (^i/e), (b) Na (^1/4) and (c) Na 
(^2/5); the layer structures correspond to the overlayer struc- 
tures in Fig. 10a, b and c, respectively. In these cumulative 
plots, the areas between successive curves show the contribu- 
tions of 5-like (red), p^+py-Wke (dark green), pAxke (light 
green) and d-\\ke (blue) states, respectively. 
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V. RESULTS AND DISCUSSION FOR W(lll)/Na 

V.l. Field and photofield emission currents from 
W(lll)/Na(lxl) 

Measurements of the TEDs from clean W(lll) at 
78 K in FE and at room temperature in FE and PFE ' 
show a strong asymmetrical peak B3 with maximum 
emission at about -0.7 eV, a broad peak A3 centered at 
about -1.4 eV^", and a rise in the experimental en- 
hancement factor of clean W(lll) just above Ef that 
suggests the presence of a peak C3. The calculated dis- 
persion plot for clean W(lll) (Fig. 12a) shows two 
bands of surface resonances of J^-like symmetry that 
originate at F at about -0.57 eV and -0.42 eV respec- 
tively. The highly asymmetrical peak B3 is attributed to 
emission from strongly lifetime broadened states B3' in 
these bands. Peak C3 is consistent with the calculated 
asymmetrical peak C3' in FE that is due to a band of in- 
termediate states of predominantly (//-like symmetry 
that originates at F at +0.18 eV, and the broad peak A3 
is consistent with a calculated weak peak A3' in FE that 
is due to several bands of surface resonances close to F . 
The energies and symmetries of the calculated peaks in 
the TEDs in FE and PFE from clean W(l 1 1) are listed in 
Table III and the energies are compared with those of the 
observed peaks. More details are given in Ref. 16. 

The bottom of the band formed by the 5-like va- 
lence states in an isolated Na layer, calculated using a 
supercell in which Na atoms occupy the same sites as in 
the W(l 1 1)/Na (1x1) supercell and the W sites are emp- 
ty, is at -1.40 eV. It disperses to higher energy along F 
M and crosses Ef at 0.75FM. At m the states in that 
band show p^+pyAike symmetry. At F the p^-Xike con- 
duction states form another band at +1.06 eV that dis- 
perses to higher energy along F M . 

Fig. 12d shows the logarithm of the calculated TED 
of the PFE current from W(l 1 1)/Na (1x1) at room tem- 
perature, plotted as a function of the initial state energy. 
The photon energy was 3.05 eV, the electric field 
strength was 0.16 V.A"\ and the work function was 
2.84 eV. The energy of the peak of the surface potential 
barrier is +1.32 eV. This accounts for the cutoff in the 
TED in PFE below about -1.7 eV. In Table III, the ener- 
gies and symmetries of the calculated emission peaks in 
FE and PFE from W(l 1 1)/Na (1x1) are listed. 

When the Na (1x1) overlayer is adsorbed on the 
clean W(lll) substrate the 5-like states of the overlayer 
hybridize with several W states near T resulting in a 
number of peaks in the K-SDOS of W(lll)/Na (1x1) 
shown in Fig. 12c. The J^-like surface resonances B3' of 
the clean substrate at -0.57 eV and -0.42 eV (Fig. 12a) 
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Fig. 12. (Color online) Dispersion plots along M F of (a) 
clean W(lll) and (b) W(lll)/Na (1x1) plotted in the surface 
Brillouin zone of the overlayer. Surface states and surface 
resonances are shown by thick curves. The predominant sym- 
metry in the surface layer is shown by the line style {s, red 
(dashed); p, green (dotted); and d, blue or grey (solid)]. Bulk 
and intermediate states are shown by thin grey dotted lines, (c) 
K-SDOS of W(lll)/Na (1x1). The successive curves in the 
cumulative plots show the contributions s-\\ke (red), Px+Pv- 
Hke (dark green), /j^-like (Hght green) and d-Wke, (blue) states, 
respectively, (d) Calculated TED in PFE for W(lll)/Na (1x1) 
with 3.05 eV photons, based on surface photoexcitation and 
plotted as a function of the initial state energy. 
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Fig. 13. (Color online) LDOS of clean W(lll) in the (a) cen- 
tral (bulk) layer and (b) surface layer, (c) LDOS in an isolated 
Na (1x1) layer. LDOS of W(lll)/Na (1x1) in the (d) W sub- 
strate and (e) Na overlayer. In these cumulative plots, the 
areas between successive curves show the contributions of s- 
like (red), p-Vke (green), fi?„+(i/.,,^-like (grey), d^z+dyAiXie 
(light blue) and d^ -like (dark blue) states, respectively. 
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Table 111. Comparison between the energies of the peaks observed in the TEDs in FE and PFE from clean 
W(l 1 1) (A3 to 03)'*''^ and the energies of the calculated peaks of clean W(l 1 1). Also listed are the energies 
and symmetries of the calculated peaks in FE from W(l 1 1)/Na (1x1) (F3' to H3'). s.r. stands for surface re- 
sonances and int. for intermediate states (states having both surface and bulk character). 



Experimental 


Calculated | 


Peak Label E-Ep (eV) 


Peak Label 


E-EpieW) 


Character 


Symmetry in 
Overlayer [Substrate] 


A3 -1.4(2)'° 


A3' 


-1.30(2) 


s.r. 


[d] 


B3 


-0.65(5) 

-0.75(5)'° 


B3' 


-0.57(2), 
-0.42(2) 


s.r. 


[dh 


C3 


+0.2(1) 


C3' 


+0.18(2) 


int. 


[d.'] 






F3' 


-1.45(2) 


s.r. 


s[d.'] 






G3' 


-0.92(2) 


s.r. 


s [d/] 






H3' 


-0.22(2) 


int. 


s [4'] 



shift to F3' at -1.48 eV and G3' at -0.98 eV (Fig. 12b) 
yielding peaks F3' and G3' in the K-SDOS. As a result 
the calculated peaks B3' of clean W(lll) that yield the 
highly-asymmetrical peak B3 observed in FE and PFE 
from clean W(lll) are suppressed. The d-like surface 
resonances A3' and intermediate states of clean W(lll) 
hybridize with the 5-like and p-like states of the overlay- 
er and yield several other peaks in the K-SDOS in the 
region from -1.70 eV to -0.98 eV. The TED in PFE in 
Fig. 12d shows a broad region of enhanced electron 
emission over that energy region due to emission from 
all these states. The unoccupied J^-like states C3' of 
clean W(lll) just above Ef shift to below Ef (labeled 
H3' in Fig. 12b) and yield peak H3' in the K-SDOS and 
in the TED in PFE. Unfortunately, there are no experi- 
mental data available in the literature on W(lll)/Na to 
compare with our calculated TEDs. 

V,2, Layer densities of states of W(lll)/Na (1x1) 

When comparing the LDOS in the central (bulk) 
layer of clean W(l 11) (Fig. 13 a) with that in the surface 
layer (Fig. 1 3b) a number of surface states and surface 
resonances are identified at about -4.0 eV, -1.0 eV, 
-0.3 eV and +1.4 eV (Fig. 13b). The LDOS of an iso- 
lated Na layer calculated using a supercell in which Na 
atoms occupy the same sites as in the W(lll)/Na (1x1) 
supercell and the W sites are empty is shown in Fig. 13c. 
The magnitude of the LDOS is constant in the range 
from -1.4 eV to + 1.0 eV, with some s-like component 
being increasingly replaced by a pt+py-like component. 
Above +1.0 eV an additional sfrong j3z-like component is 
observed. When the Na (1x1) overlayer is adsorbed on 
the clean W(lll) substrate the LDOS of the overlayer 
(Fig. 13e) is greatly modified and spreads out in energy 
to about -4.3 eV and shows several strong peaks. The 
surface resonance peaks in the LDOS of the substrate at 



-4.1eV, -1.1 eV, -0.7 eV and +1.4 eV are slightly mod- 
ified and extend into the substrate. 

VI. CONCLUSIONS 

We interpret peaks observed in the TEDs in FE and 
PFE from W(IOO) with a Na c(2x2) overlayer and from 
W(llO) with a Na (2x2) overlayer'^ by calculating the 
electronic structure and the total energy distributions of 
the FE and PFE currents from these interfaces. We show 
that Na states hybridize with surface states and surface 
resonances of clean W, shifting their energies, and re- 
sulting in various Na-induced peaks in the observed 
TEDs in FE and PFE that are consistent with experi- 
ment. We find that the strong Swanson hump Bi in the 
TEDs in FE and PFE from clean W(IOO) is suppressed 
due to Na adsorption. Instead a new peak Fi is observed 
in the TEDs in FE and PFE from W(100)/Na c(2x2) due 
to hybridization of the s-like valence states of the Na 
c(2x2) overlayer with the Jz'-like surface states of the 
Swanson hump at r'. The surface states of W(IOO) shift 
thereby by -1.2 eV. 

The different symmetry of the W(100)/Na c(2x2) 
overlayer with respect to that of the clean substrate is 
responsible for folding back the high symmetry point 

M in the SBZ of the W(IOO) substrate to r' in the SBZ 
of the overlayer, allowing states at M in the SBZ of the 
substrate to contribute to FE and PFE. p^-like conduction 
states of the Na c(2x2) overlayer hybridize with j/-like 
surface resonances of the W(IOO) substrate that have 
been folded back from M to r', shifting the surface re- 
sonances by about -0.5 eV. Bulk photoexcitation to 
these states yields peak Ji observed in the TED in FE 
from W(100)/Na c(2x2). A similar effect is observed 
when the symmetry point S in the SBZ of the W(llO) 
substrate is folded back to r' in the SBZ of the (2x2) 



overlayer. Bulk photoexcitation involving these W states 
that hybridize with the Na states is responsible for the 
strong peak K2 observed in the TED in PFE from 
W(110)/Na(5i/4). 

When a Na c(2x2) overlayer is adsorbed on 
W(IOO), the charge is redistributed among the angular 
momentum states, modifying the spatial distribution of 
charge in the vicinity of the surface. We show that the 
net effect is a strong outwardly-directed dipole layer re- 
sulting in a decrease in the work function. The work 
function we calculated for W(100)/Na c(2x2) is in good 
agreement with the experimentally observed work func- 
tion at 1 ML coverage. 

Contrary to Ishida's results**, our linear augmented 
plane wave calculations show clearly that with increased 
Na coverage from 9 = 1/6 to 9 = 2/5 the strong s-like 
peak below Ef and the strong s and j^^-like peak above 
Ep in the SDOS of W(110)/Na that are dominated by 
states at r' shift by -0.4 eV and -1.3 eV respectively. 
This energy shift is due to the increased intralayer Na 
interactions with increased coverage that broaden the Na 
states resulting in the shift of the states at r' to lower 
energy. 
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